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Abstract—In many real world scenarios, there is always a 
difference between the performances of simulated and fabricated 
structures, often due to fabrication tolerances that introduce 
statistical variations in the physical parameters, e.g., the element 
size or the periodicity of a Frequency Selective Surface (FSS). 
In this communication, we address the problem of modeling 
periodic structures with statistical variations in their geometries, 
as is typically the case with Metamaterial devices designed for 
optical wavelengths, where the difficulties in their fabrication 
almost always introduce small variations in the dimensions of 
the elements that comprise the “periodic” array. 


Index Terms—Periodic Structures, Frequency Selective Sur- 
faces (FSS), Polynomial Chaos Expansion, Latin Hypercube 
Sampling, Least-Angle Regression, Leave-One—Out Cross- 
Validation. 


I. INTRODUCTION 


Recently, there has been considerable interest in the topics of 
Metamaterials (MTMs) and Frequency Selective Surfaces (FSSs) [L]- 
[4]. These FSS elements act as spatial filters whose passbands and 
stopbands depend on the geometrical parameters of its elements. Nor- 
mally, the FSS structure is simulated during the design process and 
then fabricated to verify if indeed it has the predicted characteristics. 
It is not unusual to find that there is a considerable discrepancy 
between the simulated and measured results, even when there is only 
a minor difference between the designed and fabricated structures. 
This is especially true for Metamaterials used at optical wavelengths, 
where the difficulties in their fabrication almost always introduce 
small variations in the dimensions of the elements that comprise 
the “periodic” array. Typically, the effect of this type of variation 
in the unit cell parameters is studied by using the Monte Carlo (MC) 
methods. But if the variations are substantial, the MC method results 
in a highly computationally intensive problem. 


In this work, we use the polynomial chaos expansions (PCEs) 
in preference to the commonly used MC method, to develop a 
meta-model of the physical model of the FSS by using a polynomial 
expansion approach. In the following sections, we first describe the 
type of periodic element we are considering for this study. Next, 
we present the details of the method used to investigate the effects 
of the variations in the periodic elements. Finally, we include some 
numerical results to illustrate the application of the proposed method. 
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II. PERIODIC STRUCTURE WITH UNCERTAIN PARAMETERS 


Periodic structures are typically modeled as infinite arrays of 
scatterers, and are commonly analyzed by imposing periodic bound- 
ary conditions to a unit cell to reduce the original problem to a 
manageable size and to reduce the computational burden. In this 
study, we use cross-dipole [6J-[7] as the FSS element (see Fig. 5. 
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Fig. 1. 


The input parameters of cross-dipole FSS element 


The cross-dipole has three parameters, namely (i) the periodicity 
D+ = Dy = D; (ii) length L; and (iii) width W, as shown in Fig. 
We assume that these parameters follow normal distributions, and 
we compute the reflection coefficient for this structure for different 
values of D, L and W. 


II. THE POLYNOMIAL CHAOS EXPANSIONS 


The PCEs aim at modeling the response of a random output 
variable Y depending on input random parameters X = [11 --- Xml, 
linked by a physical model Y = M(X) [8], where the input 
parameters are supposed to be independent. Under these conditions 
we can use the PCE method to construct a meta-model of the physical 
model by using the following polynomial expansion: 


Y= Ý = aeta(X) (1) 
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Our purpose here is to substitute the physical model, which is 
computationally expensive to evaluate, by a meta-model that is much 
faster to work with. Once we have derived this meta-model, we can 
perform an MC method on this meta-model to study the propagation 
of the uncertainty of the input parameters in the physical model. 

In {ip. the Ya are multivariate orthogonal polynomials, which 
constitute a set of basis functions in the corresponding probabilistic 
space. Initially built for Gaussian input parameters, the PCE theory 
has been extended to other types of input probability density functions 
(PDFs) as the so-called generalized PCE theory [9]; for instance, 
Legendre polynomials can be used as a basis for uniform input PDFs 
spanning the range [-1,1]. 

The aq coefficients in are the unknowns that are yet to be 
computed. In the family of non-intrusive resolution methods, we 
choose the regression approach here. The principle is to truncate the 
expansion to set up a regression problem as described in , and to 
use the ordinary least-square method to compute the coefficients. The 
choice of this method is mainly motivated by its ease of resolution 
and its efficiency from the computational cost point of view [10]. 

However, since we have no information on the impact of the 
truncation on the quality of the model, this regression approach 
requires an iterative experimental procedure that uses the estimated 
accuracy of the obtained meta-model. For example, for a given design 
of experiments with a given number of simulation points, if the 
quality of the obtained meta-model is not sufficient, the design of 
the experiment should be improved with new simulation points in 
order to have a more accurate model. In this study, we will use 
the latin hypercube sampling (LHS) technique to create the 
experimental designs that we employ. To add new simulation points 
to these designs, we will use the nested LHS technique [8]. 

There are several different types of techniques that we might 
employ for the truncation of this polynomial expansion presented 
in (ip. Classically, in the regression technique, one uses the full-PC 
truncation, in which one retains the polynomials {%.,0 < |a| < p} 
for a chosen degree p. Thus, the number P of polynomials that 
are retained in the full-PC truncation depends on the degree p, and 
therefore depends on the number of points in the experimental design, 
because of conditioning problems of the information matrix arising 
in the least square estimation. For this given degree p, the number 
of polynomials kept in the full PC truncation [TO], is equal to: 


_{ M+p 
r=( M3) 2 


We observe that the number of polynomials increases with the 
dimensionality M of the input parameters and with the degree p. 
Thus, if we want a high degree p, we would need to compute a 
large number of simulation points in the experimental design. In view 
of this, we choose to employ an optimized truncation based on the 
least-angle regression (LARS) algorithm use to select iteratively 
among a candidate set of polynomials the most influencing ones onto 
the output variable that would be kept in the truncation. The so-called 
LARS truncation is described in detail in [12]. 

To assess the quality of the meta-model, the ideal case is to 
determine the generalization error. For each X, the error can be 
expressed as 


Er =E [(m% = mon)’ | (3) 


where M(X) is the physical model, M(X) is the obtained PCE 
meta-model and E is the expectation operator. But, if we only know 
M(X) for a few points (those of the LHS experimental design), we 
can at best estimate this generalization error. A well-known way to 
estimate this error is the leave-one-out cross-validation (LOOCV) 
. The principle is to take one point x out of the experimental 
design, to compute the metamodel ME? on the remaining points, 


and to calculate the prediction error at x as follows: 
AM =M (=) — 4È (x) (4) 


By computing A® for all x in the experimental design, we can 
derive the following estimate of the generalization error: 


N 
il i 2 
Errroo = y 2 A ) (5) 
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The LOO error is slightly more optimistic than the actual generaliza- 
tion error, but it yields much better error estimate than if we consider 
a simple empirical mean-squared error. From this LOO error, we can 
extract a deterministic coefficient Q? (which is the equivalent to the 
well-known R? coefficient): 


2 ErrLoo 
WY) ” 


where WY) is the estimated variance of the output Y. The closer is 
the estimate of the coefficient Q?, the better is the quality of the meta- 
model, and we will use this coefficient as a measure of the quality 
of the meta-model in this study. It has been shown in that the 
LOOCV method generally performs well in terms of generalization 
error bias. 

Here the LOOCV is used not only to estimate the accuracy of the 
generated meta-model, but also to determine the number of LARS 
polynomials kept in the truncation that provides the best possible 
accuracy. The entire procedure of PCE construction is summarized 
as follows: 


e Build an initial LHS design of experiments of a given number 
of points. 

e Perform the LARS algorithm to select iteratively the most 
influencing polynomials among a candidate set of polynomials. 

e For each possible number of LARS polynomials kept in the 
truncation, build the corresponding LARS-PC meta-model and 
assess its accuracy with the LOOCV. 

e Retain the number of LARS polynomials with the highest Q? 
value. 

e If the final accuracy of the meta-model is not satisfactory, i.e., 
if the corresponding Q? value is not sufficiently high, add 
new simulation points in the design of experiment using NLHS 
technique. 

e Restart the entire procedure from the second step until the 
obtained accuracy is satisfactory. 


IV. RESULTS 


In this study, we are interested in two different output variables, 
namely the Magnitude (Imag) and the Phase (Tphase) of the 
Reflection Coefficient. We use three input parameters (see Fig. 
that follow the normal distribution, and their properties are shown in 
Table [I] 











TABLE I 
INPUT PARAMETERS 
D L W 
Mean 45mm 30mm 6mm 
Standard Deviation | 2.25mm 1.5mm | 0.3mm 




















We have constructed a LHS design in which the simulation points 
iteratively increase from 100 by following the nested LHS technique 
[8]. For each sample of points, a best LARS-PC meta-model was 
built together with its Q? coefficient. Then if the Q? coefficient is not 
sufficiently high, we enrich the LHS design with new sample points 
until the quality of the best meta-model is sufficiently satisfactory. 
We have studied both Imag and P phase, and have obtained Q? values 
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Fig. 2. Probability Density as functions of (a) Magnitude of Reflection 
Coefficient and (b) Phase of Reflection Coefficient 


of 0.9921 and 0.9937, respectively, for Pmag and Pphase with 400 
simulation points. By taking the square root of the corresponding 
Errro9 Eb we have an estimation of the general root-mean-square 
error (RMSE) of our models. Then, the RMSE is estimated to be 
6.31 x 1075 and 3.17, respectively, for Vmag and Cpnase. With these 
two computationally easy-to-compute LARS-PC meta-models, we 
can now perform a classic MC method and deduce the distribution of 
our output variables. The resulting estimated distributions are plotted 
in Fig. 

In Fig. we can see that the Tmag values become slightly 
above 1, which is clearly non-physical. This behavior can be at- 
tributed to the sweeping window of the function used to compute 
the distribution. Indeed, since there is a high density of the value 
for 1, and theoretically none above 1, the sweeping window realizes 
a mean distribution in this area, which explains the crossing. This 
type of crossing can also be the result of the inaccuracy of the 
computed LARS-PC model. In fact, even if we have a very good 
fit on average, the error of the model can become more important in 
this region, in which there is a high density and where the output 
variable seems to have an asymptotic behavior. This is due to the 
fact that the polynomial chaos can have some difficulties to fit in an 
area where the output is locally constant. In spite of this, our model 
has very good performance in term of fitting, and we should simply 
keep in mind that values greater than 1 are non-physical. 


V. CONCLUSION 


In this communication, we have utilized a newly introduced tech- 
nique for modeling periodic structures with statistically varying 
input parameters. The proposed technique yields accurate results for 
the reflection characteristics with fewer simulation runs than those 
needed in a crude Monte Carlo simulation. 
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